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Abstract 

Using Monte Carlo simulations we study jamming and percolation processes 
upon the random sequential adsorption of dimers on binary alloys with different 
degrees of structural order. The substrates we use are equimolar mixtures obtained 
utilizing the isomorphism between an alloy and the Ising model with conserved order 
parameter. We anneal the substrates at temperature T until they reach thermal 
equilibrium, and then we quench them to freeze the state of order of the alloy at 
that temperature for the posterior adsorption trials. The annealing temperature 
is then a continuous parameter that characterizes the adsorbing surfaces, shaping 
the deposition process. As the quenched alloys undergo an order-disorder phase 
transition at the Onsager critical temperature (Tc), the jamming and percolating 
properties of the set of deposited dimers are subjected to non-trivial changes, which 
we summarize in a density-temperature phase diagram. We find that for T < 
T* = 1.22Tc the occurrence of jamming prevents the onset of percolating clusters, 
while percolation is possible for T > T*. Particular attention is focused close to 
T*, where the interplay between jamming and percolation restricts fluctuations, 
forcing exponents seemingly different from the standard percolation universality 
class. By analogy with a thermal transition, we study the onset of percolation using 
the temperature (in this case, the substrate annealing temperature) as a control 
parameter. By proposing thermal scaling Ansatzes we analyze the behavior of 
the percolation threshold and its thermally induced fluctuations. Also, the fractal 
dimension of the percolating cluster is determined. Based on these measurements 
and the excellent data collapsing, we conclude that the universality class of standard 
percolation is preserved for all temperatures. 
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I. Introduction. 

When adsorbed particles are bound to a solid surface the magnitude of the interactions 
relative to the thermal energy may lead either to reversible or irreversible attachment 
to the surface. If the particle-substrate interaction is weak enough, the particle can ex- 
plore the whole surface through surface diffusion and desorption. These mechanisms may 
eventually enable the system to reach equilibrium. The phenomenon could then be in- 
vestigated using the methods of equilibrium statistical mechanics. On the other hand, 
in the case of irreversible attachment due to strong particle-substrate interactions, the 
particles once fixed on a surface neither desorb from it nor diffuse on it. This irreversible 
adsorption or deposition of particles onto a solid substrate is a far-from equilibrium phe- 
nomenon of wide interest in physics, chemistry, biology, and in other branches of science 
and technology. Some examples embracing deposition include adhesion of colloidal par- 
ticles and proteins, separation of viruses or bacteria, adsorption of gas molecules [2, 1], 
etc. In addition to the usual case where homogeneous surfaces are used, the deposition 
can also be performed on substrates where localized adsorption takes place on particular 
sites of non-homogeneous substrates [2, 3]. The case of deposition of gas particles on a 
substrate is of great practical importance as a first step for chemical reactions in heteroge- 
neous catalysis [4] . There, the substrate properties are usually improved by alloying, and 
a large dependence of the catalytic properties on the composition and the configuration 
of the surface can be observed [5] . 

The random sequential adsorption model [1, 6] (RSA) provides an excellent description 
of the process of deposition, assuming the successive adsorption of particles within a lattice 
gas framework. Within this model, objects of finite size are randomly adsorbed on an 
initially empty rf-dimensional substrate with the restriction that they cannot overlap with 
previously deposited objects. The state of a site then changes irreversibly from empty 
to occupied. Under these conditions the system evolves with a dynamics that becomes 
essentially dominated by geometrical exclusion effects between particles. During particle 
deposition one can define different clusters looking for the sets of neighboring occupied 
sites. A particular cluster is said to be percolating if it reaches two opposite edges of the 
lattice {e.g. top and bottom). The lowest coverage at which there is a percolating cluster 
on the infinite lattice is called the percolation threshold Op. Since no desorption is allowed, 
the deposition process necessarily ends due to blocking, when no more particles fit in the 
volume; in this context we say that jamming occurs. The fraction of total space covered 
at time t by deposited particles, 0{t), reaches then a maximum value 6{t — > oo) = 9j 
called the jamming coverage. 

The saturation or jamming of a volume is an old issue that is still important today, 
linked to a wide variety of problems as relevant as car parking, occupied volume fractions 
on glasses and liquids, or packing of commercial granular goods, for which is still the 
focus of great attention [1, 7, 8]. On the other hand, percolation is one of the most 
fundamental and widely studied topics in statistical physics. The concept is applied to 
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many problems of completely different types of fields ranging from natural sciences to 
sociological phenomena. The infection of trees in an orchard [9], magnetism on diluted 
alloys [10], conductivity on complex oxides [11], and the spread of forest fires [12] are 
some popularly mentioned examples. For reviews on percolation see e.g. [9, 14, 13] . 
As in the case of thermal transitions, the percolating transition presents a non-trivial 
critical behavior, but due to purely geometrical causes. It shows scale-invariant behavior 
characterized by critical exponents with scaling relations between them. Furthermore, 
these exponents exhibit universality and do not depend on microscopic details such as the 
inclusion of next-nearest-neiglibor connectivity, or the nature of the lattice [9, 14]. 

The connection between both phenomena — ^jamming and percolation — has been at- 
tracting considerable attention [18, 19, 15, 21, 20, 22, 16, 17] and it has been shown that 
they share some similarities [21, 20, 16, 15]. In some models the deposited objects cannot 
percolate because jamming occurs before, blocking the system [23, 20]. It has also been 
shown that a continuous control parameter can be tuned to enforce the jammed system 
to go from a percolating region to a non-percolating one [24, 19, 22, 26, 25]. 

In the present paper we study the RSA of dimers on substrates composed of two- 
dimensional binary alloys [16, 17] using Monte Carlo simulations. We first obtain the 
jamming coverages and the percolation thresholds by proper extrapolation to the ther- 
modynamic limit. The temperature at which we have prepared the alloys constitutes a 
continuous variable that allows us to vary both the connectivity and the maximum cov- 
erage of the substrates. Using this temperature as a control parameter we can finely tune 
the jamming threshold in order to force its interference with the percolating transition. 
We can also generalize standard scaling concepts to a new scenario in which a thermal 
parameter intervenes in a geometrical transition [27]. All these proposals are tested by 
means of computer simulations. 

The manuscript is organized as follows: in Section II the models for the substrate and 
adsorption process are defined and the simulation method is described. In Section III re- 
sults are presented and discussed, while in Section IV the concept of thermal percolation 
is proposed, tested and discussed. Finally our conclusions are summarized in Section V. 

II. Model and simulation method. 

We study the random sequential adsorption of dimers — i.e. two identical units — on in- 
homogeneous substrates. The surface used for each deposition is a two-dimensional alloy 
annealed at temperature T and then suddenly quenched to freeze the high-temperature 
configuration. We obtained different microstates of the alloy by means of Monte Carlo 
simulations on a square lattice of side L, using periodic boundary conditions and Kawasaki 
dynamics. We took advantage of the well known isomorphism between the Ising model 
[28] and a binary alloy, namely spin-up 74-species and spin-down — > S-species, keeping 
the same density of particles pA = Pb = 1/2. The "annealing" temperature at which the 
substrate was generated is measured in units of the interaction constant (J), setting the 
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Boltzmann constant to unity. We assumed attractive interactions between species of the 
same type {i.e. J > 0, corresponding to tlie ferromagnetic Ising model). It is well known 
that this system undergoes an order-disorder transition at Tq — 2.269 in two dimensions 
[28] . Since achieving equilibrium is particularly difficult at low temperatures, we choose as 
an initial condition the microstate that minimizes the energy (two rectangular domains of 
A and 5-atoms, respectively) to save computational time. We disregarded a high number 
of configurations correlated with the initial to ensure that equilibrium was attained for 
each value of L and T. For big lattices this demanded discarding more than 10^ Monte 
Carlo steps near Tc and at the lowest temperatures. In order to perform the RSA ex- 
periments we generated and stored between 100 and 500 well equilibrated configurations 
of the alloy (depending on the substrate size) for each annealing temperature, quenching 
in this way the state of order the substrate had at this T. Subsequently, we study the 
irreversible deposition process on top of the different substrates at zero temperature — i.e. 
neglecting the diffusion of the adsorbed dimers. The only relevant temperature we will 
be refereing to throughout this work is then the one at which the adsorbing surfaces have 
been prepared. The RSA rule we assumed is the following: dimer adsorption on the alloy 
is only possible on nearest-neighbor sites with atoms of different type {AB pairs), and it 
is rejected otherwise. 

Simulations were performed by using samples of side 16 < L < 512, where distances 
are measured in lattice units. For additional details on the simulation method see reference 
[16]. 

Throughout a RSA process, the probabilities to find a percolating (jammed) cluster, 
on a finite sample of side L, can be fitted by the error function [15, 29] 

where 6 is the density of adsorbed dimers on the binary alloy, 9x{L) is its mean value 
and ax{L) is the fluctuation of that density, and x — J,P refers to the jammed or the 
percolating state, respectively. 

The binary alloy is a non-homogeneous substrate with a characteristic structure deter- 
mined by the thermal noise during the annealing period, while the adsorption of dimers 
is another random process. So, we have to deal with two correlated stochastic processes 
and the measurement of relevant physical quantities requires a careful treatment. In fact, 
if one has series of n independent samples of the substrate (i = l,..,n), it is possible 
to obtain representative values of the coverages and their fluctuations by computing the 
following averages on samples obtained at T 

1=1 

and 
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n 



i=l 



It should be remarked that the sets 6*^ and cr* , are obtained by making several adsorp- 
tion trials (~ 10^) using a single substrate and by calculating the average of the density 
and its root mean square (RMS) over these trials. It has already been shown that the 
measurement of the observables defined by means of Equations (2) and (3) captures the 
physical behavior of the adsorption process [16] (in the present case percolation and/or 
jamming). In this way, for each substrate we perform several trials searching first for 
percolation and after that for the jammed state. On other hand, the fluctuations of the 
quantity d], taken over i-indices, reflects the physical behavior of the substrate. This kind 
of study has already been performed [16], so it will not be repeated here. 

Another method for obtaining both 9"^^ and o"^ is to directly fit an error function (Equa- 
tion (1)) to the data. Although this procedure is more demanding from the computational 
point of view it will allow us to take into account the interplay between percolation and 
jamming processes explicitly in the evaluation of the relevant quantities (for example crit- 
ical exponents). 

III. Results and discussion. 



Figure 1 shows plots of the jamming coverages 6j veisu.s T obtained for samples of differ- 
ent size. These curves depend markedly on the temperature. Indeed, given the deposition 
rule that we imposed for the dimers the behavior of the saturation coverage qualitatively 
follows that of the energy of the binary alloy used as a substrate, with a steep slope at Tc 
(recall that the alloy energy is proportional to the number of broken bonds or Ai?-pairs). 
In addition to this strong variation with temperature (from near zero to above 65% of the 
entire lattice) the jamming coverage presents strong finite-size effects for T < Tq- It is 
well known [1] that the fluctuations of the jamming coverage (o"a;(L) in Equation (1) with 
the subindex x = J) scale with the lattice size according to 



where vj is the jamming exponent. A relationship similar to Equation (4) also holds for 
the fluctuations of the percolation threshold [9]. In recent papers [16, 17] we have proved 
rigorously and tested by means of Monte Carlo simulations that Equation (4) holds for a 
wide variety of RSA processes with a jamming exponent given by 
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Figure 1: Plots of the jamming coverage {9j) versus the temperature T at which the 
substrates were prepared and then quenched from. As throughout the manuscript, we 
performed the depositions considering no thermal agitation at all. The vertical dashed 
line shows the location of the ordering temperature of the underlying alloy. The coverage, 
obtained for the RSA of dimers on lattices of different size L, closely follows the behavior 
of the energy of the underlying Ising model, with a sharp slope at Tc. Finite size effects 
are more evident for T < T^, in this range of temperatures the deposition takes place 
mainly on domain walls. 

Here, D is the dimensionality of the space and df is the fractal dimension of the subset 
of active sites, i.e. sites that can allocate dimers. 

In Figure 2 the log-log plots of aj versus L show that Equation (4) actually holds 
for the whole inspected range of annealing temperatures T and lattice sizes. By fitting 
the data obtained at temperatures well below Tc {T < 2.0), we have determined uj ~ 
2/3 (see the inset of Figure 2). This result is in agreement with the idea that at very 
low temperatures the RSA process is essentially restricted to a one- dimensional interface 
between well-conformed domains of different atoms. In this case one has df — 1 and 
Equation (5) predicts i/j — 2/3 in D — 2 dimensions. 

Despite this good agreement, wc will argue that the true dimensionality of the whole 
set of sites where dimers have adsorbed is not equal to 1. When the temperature rises from 
zero it becomes increasingly probable that unlike species in the alloy start to diffuse from 
the domain wall into the bulk. They conform islands of one or more atoms surrounded 
by a sea of atoms of the other type, with their shores providing AB pairs suitable for 
the adsorption of dimers. These islands should be present for any non-zero temperature, 
increasing in number as L^. It is worth noting that these simple geometries are jammed by 
a fixed number of dimers: even though they add to the coverage they give no contribution 
to the fluctuation aj. In the low temperature range and for the sizes analyzed, we cannot 
see the large concentration of islands needed to have interference among them — or with 
the domain walls — and so the leading contribution to aj is effectively one dimensional. 
In fact, in Figure 2 we can see that above this temperature range and close enough to Tc 
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Figure 2: Log-log plots of the variance of the jamming coverage (aj) ?;ersus. L, obtained 
after deposition at different annealing temperatures T of the substrate. The inset shows 
the temperature dependence of the exponent uj obtained by fitting the curves to Equation 
(4). The redistribution of atoms in the alloy at Tc ^ 2.269 leads to a rounded step in 
Up. This exponent, in turn, allows us to evaluate an effective dimensionality for the set 
of deposited dimers. 

the interference effect starts generating a rather smooth increase of z/j when wc approach 
the critical point from below, instead of a sharp step (see the inset of Figure 2). In this 
way, df varies smoothly over effective values between one and two. 

Finally, if we look now for T > the long domain walls between A and B phases 
have disappeared (the disordering trend of the temperature has overcome the ordering 
tendency of the interactions among atoms) and one has that adsorption sites are almost 
homogeneously distributed on the sample with D = 2 and df = 2, so that Equation (5) 
yields uj — 1 in excellent agreement with the numerical results (see the inset of Figure 
2). 

In view of the previous analysis and in order to extrapolate the jamming coverage to 
the thermodynamic limit {L — > oo), we propose an Ansatz based on the assumption that 
9j{L) has two leading contributions: i) the first one corresponds to dimers adsorbed in 
the two-dimensional bulk {9j), which is independent of L; and ii) the second term {Oj""^) 
arises due to the adsorption of dimers along the interfaces between domains of different 
species, with effective dimension df which depends on L as a power law ^j"* oc L^'^f~^\ 
Then 

ej{L) = e'] + AL-^''-'^f\ (6) 

where A is a constant. Note that equation (6) resembles the scaling law generally used 
for the percolation coverage (see Equation (9) below). Figure 3 and its inset show plots of 
Oj{L) ?;ersus L~'^^'~'^f^ for various temperatures below {T < 2.20). The value of df that 
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Figure 3: Plots for the jamming coverage 9j{L) ?;ersus L'^^'^^f^ for substrates at differ- 
ent temperatures T. We took the values for the dimension df from the analysis of the 
fluctuations in the coverage aj (Figure 2). The main figure condenses some results at 
low temperature, at which df = 1, while in the inset we restrict the temperature range 
to 2.00 <T<Tc, where df departs from this value. The non-zero interception with the 
vertical axis indicates that there are dimers adding to the jamming coverage but not to 
its fluctuations. This implies that df is only an effective dimensionality for the set of 
deposited dimers. 

we used was obtained by inserting the effective exponents z/j, shown in the inset of Figure 
2, into Equation (5). The quality of the linear fits indicates that Equation (6) holds over 
the whole range considered (even for temperatures near but lower than Tc, as shown in the 
inset) according to what was discussed above, with effective dimension 1 < df < 2. The 
second term of Equation (6) vanishes for L ^ oo and the interception with the vertical 
axis provides us with an estimation of the jamming coverage in the thermodynamic limit, 
namely 9j{L — > oo) = 9j for T < Tc. We stress again that this limit would be if 
were the true dimension of the set of active sites. 

Finally, for T > Tc the jamming coverage only depends on the system size for very 
small lattices reaching a stationary value even for modest lattice sizes, as shown in Figure 
4. This fact reflects the negligible operation of lattice-size effects on the density of AB 
pairs in the bulk of the binary alloy above criticality. This finding could be anticipated 
after inspection of Figure 1 and has also been considered in the formulation of the Ansatz 
given by Equation (6) since for T >Tc one has D — df — 2. 

Percolation 

Before analyzing percolation in depth, it is worth mentioning that the maximum density of 
adsorbed dimers obtained in this inhomogeneous RSA process is very low, particularly at 
temperatures below criticality (see Figure 1). While we only have 6j ~ 0.5 for T ~ 2.80, 
this density further decreases at lower temperatures. If we consider random percolation 
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Figure 4: Linear- logarithmic plots showing finite size effects for the jamming coverage 9j 
as a function of L in a range of temperatures (T) above the ordering temperature of the 
substrate. 9j saturates for relatively small sizes, due to the two-dimensional distribution of 
the adsorbed dimers at high temperatures. Also, finite-size effects become less important 
for T much higher than the ordering temperature of the alloy (Tc). 

of dimers in the homogeneous case, the percolation threshold for the incipient percolating 
cluster is close to Op ~ 0.56 [15]. Naturally, we may note that not only the density but 
also the geometry is an important factor, and that the elongated shape of domain walls 
present in the alloy may serve to nucleate percolating clusters at low temperatures, in spite 
of the low coverage. Still, in our simulations we have found that percolation of dimers 
is not possible for T < 2.80. Of course for T = with pA = Pb = ^/'^ one has that the 
ground state of the alloy corresponds to a perfect flat interface between two domains of 
different species and on this type of substrate a trivial, one-dimensional, percolation takes 
place. However, for any finite temperature, the percolation probability of one-dimensional 
structures decreases for large samples becoming zero in the thermodynamic limit. In fact, 
when L — s> oo the probability of having a defect on the otherwise straight domain wall 
preventing the occurrence of percolation goes to 1 if T 7^ 0, which explains this result. On 
the other hand, when the temperature is increased close to T ~ 2.80 the density of dimers 
reaches a threshold that allows the onset of percolating clusters. This fact will allow us 
to draw, in the phase diagram, a percolation line starting at T ~ 2.80 that continues at 
higher temperatures. 

Let us now analyze the properties of the percolation clusters. For this purpose we will 
first test the size-scaling hypothesis for the fiuctuations of the percolation threshold and 
the extrapolation of that threshold with the system size, which are obtained in analogous 
way that we have done for 9j and its fiuctuations (see Equations (2) and (3), respectively). 
The fiuctuations (RMS) of the Op given by up scale with the system size according to [9] 
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Figure 5: Log-log plots of ap ?;ersus L at temperatures T above the critical point. While 
the fitted exponents are consistent with standard percolation for most temperatures, there 
is a significant departure on the slope of the curve measured at T = 2.80. This is the 
temperature at which the percolation line starts, with jamming and percolation occur- 
ring simultaneously in the thermodynamic limit. The inset shows 9p{L) versus L~^/'^ at 
different temperatures T above the critical point. The smallest lattice (L = 16) has not 
been included in the fit. 



ap oc L-i/'^^ (7) 

where the critical exponent vp is associated with the divergent correlation length ^ that 
behaves as 

i^\d-dp\-^-, (8) 

with vp = 4/3 for random percolation [9]. 

Figure 5 shows log-log plots of ap versus L, obtained for fixed temperatures. The 
obtained values for the exponent up are compatible with standard percolation, as it fol- 
lows from the data listed in Table I. Only the exponent measured at T = 2.80 {vp^ = 
0.80 ± 0.01) falls below the expected value, suggesting that interesting physical processes 
may take place at the point where the jamming and percolation lines meet. In order to 
extrapolate the percolation threshold to the thermodynamic limit for 9p{L) we can use 
the standard scaling approach [9] given by 

e p{L) =0p{L ^ oo) + BL-^/^p, (9) 
where S is a positive constant. The inset in Figure 5 shows plots of 9p{L) versus. L~^/^p 
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Figure 6: Log-log plots of S{L) versus L at different temperatures T above the critical 
point. We obtain the exponent Dp ~ 1.89(1) corresponding to standard percolation in all 
cases, including T — 2.80. 

obtained for T > 2.80 by taking up — 4/3, as follows from the fit of the fluctuations of 
9p{L) (see Equation (7) and Figure 5). We found that Equation (9) holds and gives a new 
confirmation of the value of the exponent vp. Furthermore, the fit allows us to extrapolate 
the percolation threshold to the thermodynamic limit for various temperatures (see also 
Table I). At T = 2.80 Equation (9) fits the data equally well for Up^ = 3/4 and Up^ = 0.80. 
Furthermore the extrapolated percolation threshold coincides, within error bars, for both 
exponents. 

On the other hand, it is well known that, at the percolation threshold, percolating 
clusters are objects with a well defined fractal dimension Dp. The number of particles of 
the spanning cluster in samples of side L {S{L)) scales as 

S{L) (X L^^. (10) 

Figure 6 shows log-log plots of the average mass of percolating clusters ?;ersus L, 
as obtained for different annealing temperatures of the substrate. In all cases (includ- 
ing T = 2.80) the results obtained by fitting the data with the aid of Equation (10) 
are in agreement with the fractal dimension of standard percolation clusters given by 
Dp = 91/48 ^ 1.896 [9] (see Table I). 

Summciry of the Results 

Figure 7 summarizes the results obtained for jamming coverages and percolation thresh- 
olds in a phase diagram. Notice that all values reported in Figure 7 correspond to extrap- 
olations to the thermodynamic limit performed with the aid of Equations (6) and (9). 
The following four regions can be observed in the phase diagram shown in Figure 7: 
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Figure 7: Phase diagram for percolation and jamming, summarizing the present results. 
We show the curves for jamming coverage 9j (open circles) and the percolation threshold 
9p (filled squares) extrapolated to the thermodynamic limit versus the annealing tem- 
perature T of the alloy. The different regions are described in more detail in the text. 
In addition to the sharp changes taking place at the ordering temperature of the sub- 
strate, there exist peculiarities at the point where the jamming and percolation lines meet 
(point JP). 
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Table 1: Critical percolation exponents measured at different temperatures as listed in the 
first column. The number in parentheses represents an estimation of the error in the last 
figure. The 2nd and 3rd columns one show the exponents obtained by fitting Equation 
(12) and Equation (10) to the simulation results, respectively. The 4th column includes 
the extrapolated values of the obtained threshold percolation by mean of Equation (9) 
using 1/ Up — 2)/ A. 



Temperature 


1/z/p 


Dp 


Op 


2.80 


0.80(1) 


1.90(1) 


0.483(1) 


3.20 


0.74(1) 


1.89(1) 


0.525(1) 


4.20 


0.73(1) 


1.90(1) 


0.546(1) 


5.00 


0.74(1) 


1.89(1) 


0.5613(5) 



Region I: This corresponds to jammed states that are inaccessible to the system. At 
low temperatures one has low coverages since the sites of the substrate suitable for dimer 
adsorption lie mainly along the interfaces between domains of different particles. By 
increasing T this scenario changes due to interdiffusion of species causing the formation 
of additional A — B pairs that implies an increase of the jammig coverages. This process 
becomes particularly relevant close to the critical temperature of the alloy (Tc = 2.269), 
so that for T > Tc one has that the jammed state is observed at higher coverages. 

Region II: Here the system has not saturated, but percolation is not observed. As 
shown in Figure 7, we have divided Region II into two parts: Region Ila, limited by the 
jamming curve and a vertical-line at T* ~ 2.80; and Region lib, above T* but below 
the percolation line. The line of solid squares above T* marks the percolation threshold. 

Region III: within this Region the coverages are high enough to observe percolating 
clusters before the system gets jammed. 

So far, along the percolation line the observed clusters belong to the universality class 
of random percolation. It seems then that the universality of the percolation process is 
not affected by the inhomogeneities of the substrate annealed at different T. This finding 
is consistent with the fact that the correlation length of the percolation process is the 
only relevant length scale. However, inspection of the phase diagram shows that there 
is a non-trivial point at the intersection of all the regions (see the "JP point" in Figure 
7). Furthermore, as we noticed before (Figure 5) the behavior of the variance of the 
percolation threshold indicates a non-standard exponent. 

Since this finding may imply a change in the universality class for percolation at this 
particular point (in what seems to be an analogy with a multicritical point in thermo- 
dynamics), we have investigated it in particular in the next section, by introducing the 
concept of thermal percolation. 
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Figure 8: Plots of the percolation probability $l(T) measured over the jamming curve 
versus T as obtained using lattices of different size L listed in the figure. ^hiT) is 
evaluated as the fraction of jammed deposition runs that have percolated for substrates 
annealed at temperature T. Solid lines correspond to fits of the numerical data to a 
Equation (11). Note that all the curves intersect at T = T* and $* = $(T*) ~ 0.93, 
implying that jamming interferes with deposition for any value of L. The inset shows 
the width of the percolating transition when we take T as the control parameter. The 
exponent we extract from the fit indicates a universality class that corresponds to standard 
percolation. 

IV. Thermal percolation. 

As it follows from the phase diagram shown in Figure 7, it is possible to cross the JP 
point from a percolating region for T > T* to a non-percolating region for T < T*, just 
by moving along the jamming curve {6j{T)) sweeping the temperature. In this way the 
percolation probability {Pl) depends on temperature, which in turn controls the prop- 
erties of the substrates. In other words one has Pl{0,T) = Pl{6j{T),T) = $£,(T), in 
contrast to the standard approach where the density is the control parameter. We will 
now explore the validity of the scaling hypothesis in this new scheme. 

Figure 8 shows typical curves of the L— dependent percolation probabilities ($l(T)) 
versus T in a range of temperatures close to T* — 2.80. It is found that curves correspond- 
ing to different sizes have an unique intersection point given by $* = ^^(T*) ~ 0.93. 

The shape of these curves, resembling so much those obtained for the percolating 
probability as a function of the occupied fraction of sites, strongly suggests testing the 
finite-size scaling approach by using the temperature as a control parameter. In order to 
do this we first fit the curves of Figure 8 by means of an error function given by (see also 
Equation (1)) 
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Figure 9: Size-dependence of the critical temperature at the percolation threshold T£ 
versus L~^/^'^ . In the fit (continuous fine) we imposed vt — 4/3. The extrapolation to 

the thermodynamic limit gives T*{L oo) — 2.80(1) 
(indicated in the graph by the arrows). 

In this way one obtains the thermal width of the transition {ctt) and the critical threshold 
T* . Of course, both quantities depend on L. We now propose that thermal fluctuations 
should scale with the size of the system in the same way as density fluctuations do, namely 
following the analogous to Equation (7). So, 

ut oc L-V-^t^ (12) 

where vt is the correlation length exponent. A log-log plot of ctt t^ersus L (see the inset 
in Figure 8) shows that the scaling form of Equation (12) holds and by fltting the data we 
obtain 1/vt = 0.74±0.01 or, equivalently, vt = 1.35±0.02. Quite remarkably, in contrast 
to the result obtained by studying percolation in the standard ensemble, the measured 
exponent for T = T* is in excellent agreement with the correlation length exponent of the 
standard percolation problem (z/p = 4/3). 

Now, the next step is to extrapolate the critical temperature — T*{L), which we have 
previously obtained using Equation (11) — to the thermodynamic limit by using an Ansatz 
analogous to Equation (9), namely 
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T*(L) = T*(L ^ oo) + CL-i/'^^, 



(13) 



where C is a constant. The obtained results are shown in Figure 9. By using Equation 
(13), with ut — 4/3 as it follows from the fit performed to the data shown in inset of 
Figure 8, we have determined a more accurate value of T* confirming that the critical 
temperature in the thermodynamic hmit is given by T*{L oo) = 2.80(1). 

Finally, the scaling laws given by Equations (12) and (13) in connection to Equation 
(11) predict the collapsing of all the curves of ^l{T) shown in Figure 8 when they are 
plotted as a function of a reduced scahng variable s = {T — T*)V'I^'^ . In fact. Figure 
10 shows a plot of the universal scaling function ^l{T) — ^(s) that results from the 
collapse of data corresponding to samples of several sizes and obtained by using the 
already determined values of both T* and ut ■ The quality of the collapse, obtained 
without any adjustable parameters, is additional evidence of the validity of the proposed 
scaling Ansatz for thermal percolation. 

In order to round out the present set of results two things remain to be explained: i) 
why we get a different exponent from the standard percolation value when the finite-size 
scaling behavior of the percolation coverage variance is analizcd close to T = 2.80; ii) 
why we recover the usual exponent when we use T as a tuning parameter for percolation 
moving along the jamming curve. 

Regarding the first issue, we will show below that the interference between jamming 
and percolation forces a different exponent when measuring at fixed T = 2.80. Indeed, 
at this temperature the threshold coverages for both phenomena are very close together 
for finite L, and they actually coincide for the infinite lattice. This implies that the 
fiuctuations in the percolation coverage are restricted by the early onset of saturation. 
In other words, in certain stochastic deposition runs either jamming occurs too early or 
percolation too late, so that the system saturates before it percolates. Then, for a fixed L, 
we are measuring a reduced value of crp, which depends on the distance \ 6j — 6p \ relative 
to (7 J. The exponent for ap that results from varying L depends then on the interplay 
of these three quantities, giving a non-trivial — and erroneous — result. Of course, this 
interference is not operative when the percolation coverage is far apart from saturation 
{Oj{L) — 9p{L) » (Tp{L)). So, in this case, we recover the normal exponents (see Table 
I)- 

In order to test these ideas we measured crp at T* in another way, trying to remove the 
constraints imposed by jamming. We measured the probability of percolation $l(^) at 
T* over 10-50 samples, and we fitted it with an error function (Equation (11)) normalized 
with the total number of runs ( 10^ — 10®). This is an important detail given that 
only a fraction of the runs has percolated (notice that the numerical data for $l(^^) — 
Figure 9 — never reach unity at any value of 9 due to the fact that ^l{9) is fixed at 

= 0.93 independent of the system size). In this way, we obtained a corrected value of 
the variance for the percolation coverage cf'p{L). After fitting this corrected values with 
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s=(T-2.80) L 



3/4 



Figure 10: Scaled plots for the percolation probability ^(s) versus the reduced variable 
s as obtained using lattices of different size L. No free parameters were fitted to obtain 
the collapse. For additional details see the text. 



Equation (7) we recovered the exponent z/ that characterizes standard percolation. We 
also tried constraining the fit of to values of 9 well below 6j, where the fraction of 

jammed depositions is negligible. Since, within error bars, we still got the same outcome 
as with the previous procedure {i.e. v ~ 4/3), we are quite confident of this result and 
the proposed explanation for the observed discrepancy. 

Let us tackle now point why do we measure = 4/3 when we move along the 
saturation curve through the percolation temperature threshold? In order to understand 
this, we will assume that when T is varied in the range of temperatures studied in Figure 
8, the connectivity changes in the underlying lattices and the adsorbed layer are not as 
important-concerning percolation-as the changes in coverage. On moving along the sat- 
uration curve (Figure 7) the coverage is fixed by the temperature T (within a deviation 
a J that is very narrow compared to up and gets narrower with increasing L [15]). We 
can then assume that T and Q are almost interchangeable or essentially linked through a 
simple functional dependence. If now we accept that Qj is linear enough as a function of T 
near T* we would then be measuring the probability of percolation in the usual way-at a 
given coverage ^j(T) -avoiding the interference between percolation and jamming. These 
arguments explain the scaling and the standard value for that we obtained from the 
thermal analysis. 

V. Conclusions. 
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Based on a numerical study of the random sequential adsorption of dimers on non- 
homogeneous binary alloys in the square lattice, we have shown that the jamming cov- 
erage (for temperatures below Tf.) and its fluctuations (over the whole range of tempera- 
tures) show the same size-scaling properties as percolation. However, the corresponding 
exponents arc different, in spite of the fact that they depend on the same dimensions 
{u — with q; = 1 for the jamming coverage and a = 2 for its fluctuations). This 

diversity is in remarkable contrast to the case of percolation, where both exponents are 
the same, being given by the divergence of the correlation length. We have also demon- 
strated that the incipient percolation cluster belongs to the universality class of standard 
percolation, as follows from the evaluated critical exponents through a finite-size scal- 
ing treatment of the numerical data. In this way inhomogeneities of the substrate are 
irrelevant for the percolation phenomena. 

In addition to these observations, an intersection between the jamming and percola- 
tion curves is found. At this particular point we observed a subtle interference between 
jamming and percolation, which seems to change one of the exponents associated with 
percolation. However, the constraint introduced by the jamming process on the perco- 
lation phenomena at this point not only does not change the universality class of the 
percolation process but, as is shown, the jamming states at different T can be used to 
characterize the critical behavior of the percolating system. In this way, we have shown 
how to use the temperature as the control parameter that governs the percolation process. 
We conclude that a generalization of the standard finite-size scaling Ansatz formulated in 
terms of the density also holds true when the control parameter is the temperature. In 
this way we are able to characterize the percolation transition at the point at which the 
two lines intersect by showing that it still belongs to the standard random-percolation 
universality class. 
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